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Abstract 

This work is a pedagogical introduction to the Lund string fragmentation model and the 
Feynman-Field hadron production model. Derivations of important formulas are worked 
out in details whenever possible. An example is given to show how to evaluate invariant 
meson production cross-sections using the Monte Carlo integration technique. 
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1 Introduction 



The topics discussed in this work are found in original sources hsted as References [1, 2, 3]. 
The materials shown here are parts of the appendices in my Ph. D. dissertation [4]. When 
I was going through these materials the first time, 1 found some typos and occasionally 
inconsistent notations in the original works. There were also some elementary questions 
that one would normally ask, which are not explained in the texts. This paper is essentially 
a record of my notes and reflections when I was learning these topics. The Lund model and 
the Feynman-Field model are chosen for this discussion because they form the theoretical 
foundation of the latter part of my dissertation when calculating the invariant hadron 
production cross sections in the soft and hard Pt regions respectively. The core of the 
latter part of my Ph. D. dissertation includes (1) a non-perturbative calculation of the 
invariant meson production cross sections in the soft pr region based on the Lund model, 
(2) a perturbative calculation of the invariant production cross sections in the hard pr 
region based on the Feynman-Field model, and (3) the parameterizations of the invariant 
pion and kaon production cross sections in both the soft and hard pr regions. The new 
results listed above will be published as a separate paper and not to be discussed in this 
work. 



2 The Area Law of the Lund Model 

The Lund model is discussed extensively in Reference[l] along with many other interesting 
topics. We shall focus only on the Area law in this section. 

2.1 Kinematics 

The quark-antiquark pair of a meson is massless in the string model. In this case, the 
concept of the mass of a meson is associated with the mass of the string field and not the 
quarks. Massless quarks move at the speed of light. It is not a surprising result given 
the fact that string theory predicts that the ends of open strings move at the speed of 
light either by imposing a Neumann boundary condition for open strings [5] or by solving 
the classical equations of motion [6]. In the highly rclativistic problems, the light-cone 
coordinates are the natural choice. In the (1-1-1) case, Lorentz transformations can be 
written as 




cosh y — sinh y 
— sinh y cosh y 



t' \ _ ( 7 —^v \ I t \ _ I cosh y — sinh y 
^' j \ 7 / \^/~\~ ^i^^ y cosh y 





(1) 

(2) 
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where 7 is the Lorentz factor 
These equations imply 
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7 = cosh|/p, 
= sinh Hp, 
V = tanhyp. 

The energy and momentum of a particle can now be written as 

E = ■ym = m cosh 

p — ^mv — m sinh y^. 
The boosted energy Ei, and the boosted 1-momentum are given as 

Eb = -f{E-vp) 

— m (cosh Up cosh y — sinh yp sinh y) 

= m cosh(|/p - y), 
Ph = 7(p-w^) 

= m (sinh yp cosh |/ — cosh yp sinh |/) 

= m s\nh{yp- y). 

The relativistic velocity addition formula is simple in light-cone coordinates [1] , 

v' = tanhy' = ^ — — = tanh(7/ - yb), 

1 — VbV 

illustrating the additivity of rapidity, 

y^y' + Vb- 

This simplication motivates the definition of rapidity 

1 /l + v\ 1 fE + p\ 



2 \l-vj 2 \E-pJ 
Momenta along the hght-cone can be defined as [1] 

p^ = E + p — m cosh yp + m sinh yp — m e^^ , 
P- — E — p — m cosh yp — m sinh yp — m e~^^. 



With these definitions, boosts are simphfied along the hght-cone [1], 

p'^ = m e±(2/P-2/) = p^e^^. (16) 
Light-cone coordinates can also be defined in configuration space as 

t = — cosh^/g, (17) 

TTl 

X = — sinh7/g, (18) 

K 

where k is the string constant and is used here to give t and x the correct dimension. 
A new subscript is adopted for Uq in configuration space to distinguish Up in momentum 
space. These definitions are consistent with the requirement that v = dx/dt — tanhy^, as 
in Eq. (6). Similar results are obtained as in the momentum space case, 

x^ — t -\- X — — cosh yg^ smh yq — — e*''^ , (19) 

K K K 

m ^ m ^ m _., , ^. 

X- — t — X — — cosh yn smh yn — — e (20) 

K K, K 

= - e±(^«-^) = x±e^y. (21) 

K 

2.2 Deriving the Area Law 

The focus of this paper is primarily mesons. The force field between a quark- ant iquark 
pair is presumed to be constant and confined to a flux-tube. The equation of motion of 
any one member of the quark- antiquark pair acted upon by a constant force is [1] 

where k is the string constant. The solution of the equation is [1] 

p{t) =po- nt = n{to-t). (23) 
Prom E'^ + m^, EdE dp or 



p _dE 
E~~dp 



(24) 



is obtained. There is also ^ ^ 

p = ■ymv = ■ym — = E —. (25) 

CLv CLL 
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Together, Eqs. (24) and (25) yield [1] 



dx p dE 

M^E = ^^ ^^'^ 

such that [1] 

dE dE dp dt dp 

dx d,p d,t dx dt 

Its solution gives the expected QCD flux tube result [1] 

E = k{xo-x). (28) 

Combining Eqs. (23) and (28) and m — 0, the relativistic equation of state is [1] 

m^^E'^-p'^^ k[{xo - xf - {to - tf] = 0. (29) 

Let Xq = to = 0, the equation of motion |a;| = t is obtained. The maximum kinetic energy 
at a; = puts an upper bound on the displacement of the particle such that nXmax = ^ ■ 
The path of the particle is illustrated in Figure 5. The zig-zag motion of the particle is 
the reason for the name "yoyo state." The quark-antiquark pair of a meson are assumed 
to be massless. The mass square of the meson is proportional to the area of the rectangle 
defined by the trajectories of the quark-antiquark pair. The period r of the yoyo motion 
is [1] 

r = , (30) 

K 

where Eq is the maximum kinetic energy of each quark. In a boosted Lorentz frame, the 
period is transformed as [1] 

r' = r cosh I/. (31) 

The breakup of a string occurs along a curve of constant proper time such that the process 
is Lorentz invariant. The Lund model assumes that the produced mesons are ranked, 
meaning that the production of the n-th rank meson depends on the existence of rank 
n— 1, n— 2, . . . , 1 mesons. A vertex y is a breakup point and the location «;(x+, xJ) where 
a quark-antiquark pair is produced. The breakup of a string is represented in Figure 6. 
The produced particles closer to the edges are the faster moving ones, corresponding to 
larger rapidities. 

Let p±o and p±j be momentum of the parent quark and the j-th of the n produced 
quarks moving along the x± light-cone coordinate respectively. Then [1] 

n 
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and the momentum fraction at Vj is defined as [1] 

= (33) 

P±Q 

In order to simplify notations, Zj — z+j unless specified otherwise. The invariant interval 
is 

ds'^ = dt'^ - dx^ = dt^{l- -y^), (34) 

giving 

ds = dtVl^ ^ — ^dr, (35) 
7 

where r is the proper time. The proper time is also x^x_ = {t — x){t + x) = — x"^ — 
and can be used as dynamic variable such that [1] 

T^K^x+X-. (36) 

Since V is proportional to the proper time square r^, its Lorentz invariant property makes 
it a good kinematic variable in the light-cone coordinates. Let W±i and W±2 be the 
kinetic energies along x± at vertices 1 and 2 respectively. The following identities can 
be easily proven geometrically by calculating the areas of the rectangles Vi = Ai + A3, 
V2 — A2-\- A2, and w? as shown in Figure 7 [1] , 

Ti = {l-z^)W-2W+i, (37) 
V2 = {l-z+)W_2W+y, (38) 
= Z-Z+W-2W+1, (39) 

where m is the mass of the produced meson. Eqs. (37) and (38) can be expressed with 
the help of Eq. (39) as [1] 

(40) 
(41) 





m\l- 




Zj^Z- 




T2 = 


m\l- 




z^z^ 




;ives [1] 






9ri 


2 

m 




dz- 


z^zt 




dT2 


2 

m 




dz^ 


z\z_ 





(42) 
(43) 



which will be used later in the next section. 
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Let H{ri) be a probability distribution such that i^^(ri) dTi dyi is the probabiUty of 
a quark-antiquark pair being produced at the spacetime position Vi. From now on, the 
symbol Vn also represents the breakup event that leads to the creation of the n-th quark- 
antiquark pair along a surface of constant r. Let f{z+)dz^ be the transition probability 
of obtaining V2 given that Vi occurs. The transition probability of V2 via Vi is then 
H{ri) dVi dyi f{z^)dz^. Similarly the probability of Vi via V2 is H(r2) dr2 dy2 f{zJ)dz-. 
It is reasonable to assume that the probability of Vi via V2 is equal to that of V2 via Vi 
such that [1] 

H{T^) dFi dyi f{z+)dz+ = H{T2) dV2 dy2 f{z^)dz^. (44) 



The Jacobian J in 



IS 



J 



dTi dz+ — 


JdV^, 












dz- 




dz+ 






dT2 


9ir 



Z- ' 



(45) 

(46) 



which can be easily computed by using Eqs. (42), (43) and dri/dV2 — dz+/dz- — 0. 
Eq. (45) can now be re-expressed as [1] 



dz^ 



dr. 



dz- 



z^ 



Since rapidity is additive (i e. 2/2 = yi + ^v)-, 

dyi = dy2. 

With Eqs. (47) and (48), Eq. (44) is simphfied as [1] 

H{Vi)z+f{z+)^H{T2)z^f{z^). 
New definitions are now made to facilitate the solution of the equation [1], 

Mr) - iogi/(r), 

g{z) = log(^/(^)). 
The new definitions transform Eq. (49) as [1] 

h{V^)+g{z+)^h{V2) + g{z^). 

Notice that [1] 



dz+dz- dz+dz- 



= 0. 



(47) 



(48) 



(49) 



(50) 
(51) 



(52) 
(53) 
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Differentiate Eq. (52) with respect to 2;+ and z^^ eliminate terms with d'^gjdzj^dz^ using 
Eq. (53) and cancel out factors of dV /dz on the both sides of the equation to obtain [1] 

dh{V,) , ^ d^hjV,) _ dh{T,) , ^ d^hjV,) 

^fT+^^^fr-^fr+^^^fT' ^^^^ 

or equivalently [1] 

where 6 is a constant. The solution is [1] 

/i(r) = -6r + alnr + lnC, (56) 

where C is a constant of integration. It yields the distribution [1] 

//(r) ^e'*^^) = Cr«e-^^. (57) 

Substituting Eqs. (40) and (41) into Eq. (52) gives [1] 

. . hm? , , 1 - ^+ , ^ 
gi2\z+) H h Oi in a2 m h m Gi 

z+ ^+ ^+ 

})Tn^ Tfl^ 1 — z 

^g2l{z-)^ ha2ln oiln ^ + lnC2, (58) 

Z- Z- Z- 

where 5'i2(^+) is g{z+) of Vi transitioning to V2 and g2i{z-) is g{z-) of V2 transitioning to 
Vi. If a = ai = a2, the transition probability distribution is [1] 

1 6m^ 

/(^,-) = ^-(l-^.re"^^- , (59) 
^3 

where is a constant of integration. Otherwise, the transition probability from to V/3 
is [1] 

fa/3(zj)^N^p-z]'^ [ e , (60) 

where N^t^ is a constant of integration specific to the vertices and Vg. Let z^j be the 
j-th rank momentum fraction scaled with respect to p+o, then [1] 

zoi = zx, (61) 
^02 = 22(1-^1)- (62) 

The probabihty of two dependent events is the product of the probabihties of the two 
individual events. The existence of the rank-2 hadron depends on that of the rank-1 
hadron according to the Lund Model so that joint probability of their mutual existence 
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is the product of the two individual probabihties of Vi and V2. With Eqs. (60-62), the 
combined distribution of the 1st and 2nd rank hadrons is [1] 

f{zi)dzif{z2)dz2 = f{zoi)dzoif f ) ^ 

V J- — ^01 / J- — ^01 

(1-^01) 1-:; e -1 -2 



zoi Z02 V 1 — 2:01 

= ^^^^(l-^oi-^02)'^e-''(^^+^^). (63) 
^01 ^02 

The geometical identity Aj — bm^/zj is used in the last step. Generalizing the product 
of two vertices to that of n vertices, the differential probability for the production of n 
particles is easily seen as [1] 

dP{l, ...,n)^{l-zrfl e-"^^ , (64) 

where z — Z^"=i Zoj. Let poj — Zojp+o and d^p — dp+ dp^. Use the identity [1] 

j dCdB5{BC -D)^^ (65) 

to obtain [1] 

d Poi Po2^"^(Poi - ) ^"^(Poi -m). (66) 



^01 ^02 

With Eq. (66), Eq. (64) can be rewritten as [1] 

n 

dP(poi, . . . ,Pon) = (1 - n ^^Vi - e"'""^- (67) 

It can be shown from simple geometry in Figure 8 that the kinetic energies of the quarks 
along the ± light-cones are [1] 

= ^Po+, (68) 

^ 779 ^ 

and the total kinetic energy square at Vn is [1] 

s = W^W_ = jz'^- (70) 
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The total differential probability of n-particle production is [1] 
dP{z,s;pou ■ ■ ■ ,Pon) 



dz5 \z - Y^zqA ds5\s-Y^ dP{poi, . . . ,pon) 

V j=l J \ 3=1 

— dsd Is dP{poi,...,pon), (71) 



where uj = po+j/Wn+- Use the trick [1] 



j=i ^ I V j=i 



j=l ) \ j=\'^j^n+^ 

where Prest = iWn+i W^n-)) to re-organize Eq. (71) as [1] 

dz / " 

dP{z,S]PQi,. . . ,PQn) = ds {1- Z)" e'''^6 \prest-J2P0j 

n 

xUN d^poj 5+{pI^ - m") e-^^--, (73) 

3=1 



n 



with the definition [1] 

Atotal = J2^3 ='^ + ^rest- (74) 
3=1 

The claim is that Afotai is Lorentz invariant. It is called the "Area Law" [7]. Eq. (73) can 
be separated in the external and internal parts as in [1] 

dz 

dPea^t = ds-{l-zre-'"', (75) 
z 

dPint = flNd^poj5+(pl^-rn^)5lprest-j2Po3]e-''^''''*. (76) 

3=1 \ 3=1 J 

The external part contains kinematic variables s, z and F. The internal part contains 
dynamic variables poj- Eqs. (75) and (76) are the final results. 
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3 The Parton Model 



The relevant formulas of the parton model in reference [2] are derived below. Several 
mistakes found in reference [2] have been corrected here. 

3.1 External Cross Section Formulas 

The Feynman diagram of hadrons A and B scattering into partons c and d, A + B c+d, 
is illustrated in Figure 9. The external Mandelstam variables are defined as 

s = (Pa + Pb)^ = 2Pa-Pb, (77) 

t = {p.-PaY = -2p,-Pa = ~xtZ, (78) 

u = (p^-Pb)^ ^-2pc-Pb^-^xt^, (79) 



where 



and 



2 % 



XT ^ (80) 

Q 

7; = tan— . (81) 



2 

Eq. (78) is proved as follow: In the high energy limit, rrii = 0, Vi, such that Ei = |p 



Ea = Eb = \Pa\ = \Pb\ = (82) 
From Figure 9, it is obvious that 

Pt = I Pel sin 6c. (83) 



-2Pc-Pa 

-2{E,Ea-\Pc\\Pa\cos9,) 

-2 (IPcIIP^I - |Pc||Pa| cos 6^0) 
, 1 — cos9r 

-2pt|Pa 



-2 



sin 9c 

(y^\ l-(l-2sin^(ge/2)) 
\2 J 2sm{9j2)cos(9j2) 



(84) 
(85) 
(86) 

(87) 



(89) 
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Similarly, 

u = -2p^-Pb (90) 

= -2(£;e£;B- |Pe||PB|cOs(l80°-^c)) (91) 

= -2(£;e£;B + |pc||PB|cOS^c) (92) 

= -'^^T^. (93) 

In the massless limit, the sum rule of external Mandelstam variables is 

s + t + u = 0. (94) 

A fraction, Xi, of the momentum of the incoming hadron Pj makes up the momentum of 
parton pj. More specifically, we have 

Pa = XaPA, (95) 

Pb = XbPs- (96) 



These imply that 

|P.| = (97) 

|P.I = (98) 

The internal Mandelstam variables, s, i and -u, of the partons can now be written in terms 
of the external Mandelstam variables, s, t and u. 

S = (Pa + Pbf = '2-Pa ■ Pb = XaXbS, (99) 
i = (Pc - Paf = -'2pc ■ Pa = Xat, (100) 
U = (pc - Pbf = -2pc • Pb = XbU. (101) 

The conservation of energy-momentum is given as 

Pa+Pb=Pc+Pd- (102) 

Prom the expressions of the internal and external Mandelstam variables, it is easily seen 
that 

o tU tU 

Pt-- = -- (103) 
In the massless hmit, the sum rule of the internal Mendalstam variables is 

5 + ^ + ^ = 0, (104) 
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or equivalently, 

From Eq. (105), we have 



XaXbS + Xat + XbU = 



—XbU —XbU/s XbXi 



XbS + t Xb + t/s Xb — X2^ 



where 



Xi 
X2 



U 1 1 

* 1 -r 

s 2 



Similarly, 



Xb = 



S XaX2 



XaS + U Xa + U/S Xa - Xi 

The expression t — {pc — PaY — {Pd ~ PbY be considered separately as 

i = iPc- Paf = -'^■Pc- Pa = -lj^XTXa%, 

i = (Pd-Pbf ^ -'^Pd-Pb^ —^TXb^, 

^ -id 



where 



Eqs. (110) and (111) imply that 



Td = tan 

a 2 



•^a — ■^b ■ 
J-d 



(105) 
(106) 

(107) 

(108) 

(109) 

(110) 
(111) 

(112) 
(113) 



The internal Mandelstam variables can be re-expressed in terms of the external angles as 
follow. 



|Pa-Pb| = |pc| COS^c + |Pd| COS^d 

Pt /) , a 



sin 9, 



sin 6*^ 



(114) 
(115) 



which implies 



\xa Pa — Xb Pb| 



PT(cot Oc + COtOd), 



(116) 
(117) 
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or equivalently 



XT (cot 9c + cot 9d) 

/I -tan2(^,/2) 



XT 



XT 



2tan(^e/2) 



+ 



2tan(^d/2) 



27; 



+ 



2X, 



— - Xt 

2 



Eqs. (113) and (121) together give 



Xb 



1 (1 1 

^Xt {% + %). 



(118) 

(119) 

(120) 
(121) 

(122) 
(123) 



Finally, the internal Mendalstam variables can be re-expressed in terms of the external 
angles as 



S — ~ Xn 



u 



% 



(124) 
(125) 
(126) 



Eqs. (124)-(126) will not be used directly in this work. However it is interesting to know 
that the internal Mendalstam variables can be measured in terms of the external angles. 
Another interesting fact worth mentioning is that Eqs. (106) and (109) are still true even 
if Tc is replaced by 7^ in Eqs. (107) and (108). This statement can be checked easily by 
substituting Eq. (113) into Eqs. (106) and (109) and then solving for Xa and Xb- The 
interchangability of % and 7^ will become useful later in some of the derivations below. 
The differential cross section of parton production for the reaction A + S— >c-|-dis 



da 



da = fA/a{Xa, Q^)dXa fB/b{Xb, Q'^)dXb dt, 

dt 



(127) 



where fA/ai^a,Q^) and fB/b{xb,Q'^) are the renormalized parton distribution functions 
(PDF's) of parton a in hadron A and of parton b in hadron B respectively. Although 
has a simple interpretation in QED, 



-E'Esm- 
2 2 



(128) 
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(where p and p' are the incoming and OTitgoing lepton momenta and 6 is the scattering 
angle), there is no clear understanding of what is in hadron collisions. Some good 
guesses are 

= iPT, (129) 

Although the interest of this work is not parton production per se, it serves as an intuitive 
starting point to introduce hadron production later. 

The substitution of Eqs. (107) and (108) into Eq. (106) gives 

1 11 

XaXb - 2 XaXrX: = ^ X^Xt (131) 

Differentiating in Eq. (131) with respect to 7^ and xt while keeping x^ fixed gives 

dX], XaXxTc - XbXT 



8% ^Xa% - XT 

dxb _ Xb + XgT^ 

dxT '^Xa'Tc — Xt 

Repeating the same process with t — —^xtXoXc yields 



(132) 
(133) 



di S 

^ = -^XTXg. (134) 

^ -^Xa%. (135) 



dxT 2 

The Jacobian now can be computed with Eqs. (132)-(135) as 

d{Xb, i) Xa Xb Xt S Xa Xb Xt s 



or 



d{%,XT) 2Xa%-XT 2{Xa-Xi)Tc 



d{xb,t) d{xb,t) d% 



d{ec,XT) d{%,XT) dOc 

Xa Xb Xt S 



2{Xa - Xi)Tc \ 
S XaXb Xt 



2 Xa — xi sin 9c 



(136) 

(137) 

(138) 
(139) 
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Given that E = ttit cosh y and Pz = mx sinh y, dpz/ dy = E. In the cyhndrical coordinates, 



dp^ = ndpzdp^ (140) 
= nEdydp^. (141) 



With Px = dxT, 
And with 
there is also 



dpr = - xt dxT- (142) 

1 l + COS^c (^.^. 



sin 6^ 

Substitute Eqs. (142) and (144) into Eq. (141) to get 



^ = ^dydpl (145) 
= "^^de^dxT. (146) 

Combine the Jacobian in Eq. (139) with Eq. (146) and obtain 

dx.di^-^^^^. (147) 

TT X(x — X\ E 

With Eq. (147), the Lorentz invariant differential parton production cross section from 
Eq. (127) can be written as 

^% = I f^"^^- ^"/'t-- f ■ 

The inclusive cross section for A+B c+X can be obtained from Eq. (148) by integrating 
over 6fi, or equivalently by integrating over Xa as in dxa = {xa — Xi) dy^. The relationship 
between dyd and dxa can be derived as follow. Substitute Eqs. (113), (107) and (108) into 
Eq. (109) and obtain 

Xad XtTo) 

Xa%%^ ^\ ^ . 149 
Xa can be isolated from Eq. (149) as 

11 , , 

Xa - 7,Xt^ ^Xi. (150) 
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Differentiating Eq. (150) witli xt and xi lieid fixed yields 



dxa = ]-XT^d% (151) 

11 , , 

= -XTTj^dya (152) 

= -XT^^dya (153) 

= ^c^y, (154) 

Xb 

= (xa - xi) dya. (155) 
Tlie inclusive cross section for A + S — > c + X is 



, 3 - - / C?2;„ /A/a(2;a, <3^) /b/6(2;6, Q^) ^"^'^ (156) 

where 

xf" = -^ = -^^/%. (157) 

The lower limit of the integral in Eq. (156) is not zero is that Xa, as seen in Eq. (106). 
Xa < for Xb < X2, which is not allowed. At which is also impossible. 

Hence the value of Xa is limited by a range of Xb- It can be seen in Figure 10 that Xa is 
at its minimum when Xb = 1 and the maximum value of is 1 . 

Once the formalism of parton production is complete, the derivation of the inclusive 
hadron production cross section can be essentially the same. The external Mandelstam 
variables are similarly defined, 

s = {Pa + PbT = '2Pa-Pb: (158) 
t = {Ph-PAy = -2ph-PA = ~XT%, (159) 

u = {Ph-PBf^-2ph-PB^~XT^. (160) 
The internal Mandelstam variables are slightly modified as 

S = {Pa+Pbf ^ XaXbS, (161) 
i = {pc-Pa)^ = -2pc-Pa = -2—-Pa=—: (162) 

u = (Pc-PbY ^ -2pc-Pb^ -2— -Pb^—, (163) 
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where 



Pa — ^aPA: 
Ph = ZcPc- 

Eq. (104) combined with Eqs. (161)-(166) gives 

X2 Xi 

I" ) 
Xb Xd 



with 



u 1 1 

s 2 % 
t 1 ^ 

X2 = = -XTTh, 

s 2 
7/j = tan-^. 



Differentiate Eq. (167) to obtain 



dzc 




dOcrn 




dzc 


1 / 1 


dxT 


2 \Xa% ^ 



H sec 

Xb J 

% 

Xb 



dt^ _ 1 XTXg 2 ^cm 
de ~ 4 ^^"^ 2 ' 

di ^ 

2 

The Jacobian can be computed from Eqs. (171)-(174), 

d{zc,i) xts 



d{Ocm,XT) 2ZcSmecm' 

As before, 

2 sin Ocm dp^ 

dOcrndXT = 

TTS Xt ti 

such that ^ 

ci^^c cii = 
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The differential liadron production cross section is 

do- = fA/a{Xa,Q^)dXafB/b{Xb,Q'^)dXb^diD'^{Zc,Q'^)dZc. (178) 

The major difference between the parton and hadron cross section formulas is the inclusion 
of the fragmentation function (FF), D'^{zc,Q^)- The inclusive hadron production cross 
section is 

E—^- dXa . dXbfA/a{Xa,Q')fB/b{Xb,Q')D^^{z,,Q^)--^, (179) 



with 



i- — X2 
Xa X\ 



3.2 Internal Cross Section Formulas 

The internal cross sections, 



^{ab^cd;s,i) ^ """"ly^ \M(ab^cd)\\ (182) 
dt 

are derived by adding the leading order Feynman diagrams [13] together, with 

4 s"^ -\- 

\M{mj ^ mj)]' = IMiqi-qj^qiqj)]'^-^^, (183) 
^ + (184) 

|A<to-«W.)r = g^— j-3(— j, (187) 

m9 9 -Ut MS St\ 

IM9S-Ss)r = 5 , (189) 
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and 

^^0{Q') = . J-^^ .Q,y (190) 

(33 - 2nf) log (^J 



1 306 - 38n/ .^^m , f Q^^ 

«L0(<3 ) log log ' 



(191) 



47r 33 -2n/ ""'^ ' ° ° \A\ 
where n/ is the number of flavors and Ajjg pa 0.2 GeV^ by convention. 

3.3 Fragmentation Functions 

The fragmentation function, D'^[z,Q'^), registers the number of hadrons of type h with 
energy fraction, 

2Eh 



(192) 



^" Q ' 

per dz generated by the initial quark q. Energy conservation gives 

f'zD'^{z,Q')dz^l. (193) 
h ""^ 

It is assumed that scaling occurs, i. e. 

D'^iz,Q')^D'^{z). (194) 
The multiplicity of h emerging from q is 

/' D'^{z)dz, (195) 

where Zmin — Sm^/Q. The fragmentation function cannot be calculated exactly by QCD. 
Feynman and Field have a semi-analytical parameterization [2], which is now being de- 
scribed below. Let /(?]) drj be the probability that the flrst hierarchy (rank 1) meson 
leaves fractional momentum t] to the remaining cascade. The function /(?]) is normalized. 



f' f{7])d7]=l. (196) 

Jo 



Secondly let F[z) dz be the probability of finding a meson (independent of rank) with 
fractional momentum z within dz in a jet. Feynman and Field assumed a simple stochastic 
ansatz that probability distribution is recursive. This way F{z) satisfies the integral 
equation 

Fiz) = /(I -z)+ fiv) F{z/v). (197) 

Jz r) 
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The first term gives the probability that the produced meson is rank 1. The second term 
calculates the probability of the production of a higher rank meson recursively. Data of 
deep inelastic scattering experiments suggest the simple form 

zF{z) = f{l-z) = {n+l){l-z)''. (198) 

The power n — 2 gives a qualitative description of experimental data. 

Next the flavor dependent part of the fragmentation function is considered. Let (3i be 
the probability of the qiq, pair. The normalization condition imposes that 



EA = 1- (199) 

i=l 

The isospin symmetry implies that 

Pu^Pd^ p. (200) 

Furthermore, data indicate that /?« ~ | /5« and that (5c and are small. For a quark of 
flavor the mean number of meson states of flavor ah at z is, in analogy of Eq. (197), 
given by 

Pf{z) = 5„,A/(1 -z)+ f'^ f{rj) (3, Pf{z/ri). (201) 
The mean number of meson states averaged over all quarks is 

Pi>{^)-Y.(^cPf{zl (202) 

c 

so that 

Pi\^{z) = -z)+ £ ^ fiv) Pt^{zh). (203) 

Comparing with Eq. (197) yields 

Pi\^{z)^(3a(3,F{z). (204) 
Comparing Eqs. (201), (203) and (204) yields 

Pf{z) = 5aM{l -z)+ MbF{z), (205) 

where 

F{z)^F{z)-f{l-z). (206) 
Eqs. (198) and (206) together give 

zF{z)^{n+l){l-z)''+K (207) 
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The mean number of meson state, (z), is related to the fragmentation function as in 

^J(^) = Er;^n"'(^), (208) 

a,b 

where F^^ is the probabihty of h containing ab. For example, F^^ = 1 and FJ^^ = F^^ = |. 
Substitute Eq. (205) into Eq. (208) to obtain 

D'^{z) = A'j{l-z) + B'^F{z), (209) 

where 

K = EKl^b^ (210) 

b 

= Y.Mr (211) 

As a example, given that /3u — Pd — 2/3s = /3 and Pc — Pb — Pt — 0, it is easily shown 
that /3 = 0.4, and 

(212) 

(213) 

(214) 
(215) 





2^' 


A-: 


= ift 

2^' 


a: 


= 0, 


B'" 


= p'. 



3.4 Implementing the Fey nman- Field Model 

Feynman and Field calculated the invariant production cross section formula by incorpo- 
rating the structure functions (or pardon distribution functions), fA/ai^a-iQ"^)-! obtained 
from DIS experiments and the fragmentation functions, D^{z, Q"^), derived from a combi- 
nation of stochastic arguments and parameterizations of data. The cross section formula 

fl^n 1 /■! /■! 1 rjfr 

E-r^ = -J2 dXa . dx, fA/a{Xa: Q') fB/b{x,, Q^) D^,{z,, Q^) - , (216) 

dp-^ TT ~^ Jx^'^ Jx]^'"^ Zc dt 



with 



C'" = T^, (217) 

L — X2 

^min ^ XaX2 ^ ^218) 
X(j Xi 
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is derived explicitly in Appendix 3. The invariant cross section is of interest because it can 
be used to calculate a number of physical quantities. For instance, the particle number 
Mh of hadron h can be calculated as [8] 

Mp{l + K,,) rf'V , , 



- j Vir) E ^ Wip, r) ^ d\ (219) 

where Mp is the protons incident on the target; W{p, r) is the probability of recording a 
pion produced in the target with coordinate r and momentum p, V{r) d^r is the probability 
for the incident proton to interact in the element d^r in the target; K^, is a coefficient 
which takes into account the probability of pion production as a result of two or more 
interactions in the target. Blattnig et al. [11] calculated the spectral distribution and the 
total cross section from the invariant cross section as in 

^ = 27rp / dOE^ sin (220) 
dE ^ Jo dp^ ' ^ ^ 



and 



(7 = 27r/ de dpE — -^ , 221 

^0 Jpmin dp-^ /p2_^^2 

where ruh is the mass of the produced hadron. These are just some examples to show 
why Ed^a/dp^ is an interesting quantity to calculate. 

Monte Carlo integration package VEGAS is used to calculate Eq. (216). VEGAS uses a 
mixed strategy combining importance and stratified sampling. The details of this method 
can be found in references [3, 9]. Inputs used in VEGAS are given as follow. 



ndim 


= 2, 


(222) 


ncall 


= 5000, 


(223) 


itmx 


= 5, 


(224) 


nprn 


= 0, 


(225) 


xl[l] 


= regn[l] = 

I - X2 


(226) 


xl[2] 


roi ^2 
= regn[2\ = 

1 — Xi 


(227) 


xu[l\ 


= regn[3] = 1.0, 


(228) 


xu[2] 


= regn[4] = 1.0. 


(229) 



The dimension of the integration is ndim. The number of random calls in the Monte 
Carlo integration is ncall. The maximum number of iterations is itmx. The fiag nprn 
controls the initial conditions of the grids. In this case, nprn = signals a cold start. The 
random sampling is confined to an n-dimensional rectangular box. In Lepage's original 
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paper, the coordinates of the "lower left corner" are labelled as xl[i]. and those of the 
"upper right corner'' are xu[i]. For a 2-dimensional integral, i G {1,2}. Numerical 
Recipes adapted Lepage's code [3] and puts xl[i] and xu[i] in an array labelled regn\j] 
where j E {1, • • • , 2 * ndim}. The integrand of Eq. (216) is implemented as the function 
f xn and is expressed in pseudo-code as 

if (x_b < x_a*x2/(x_a-xl) ) return 0; 

else (sum all the terms in the integrand) ; 

The inclusive cross section of pp — > tt+X is chosen to illustrate the Monte Carlo integration 
method in this section. The most accurate parameterization of the parton distributions of 
proton is given by the CTEQ6 package [10]. The QCD running coupling constant, as{Q'^), 
used in the code is the renormalized coupling constant described in Appendix 3. A typical 
value of A = 0.4 GeV is used inside as{Q^)- The internal scattering cross sections of the 
reactions qiq^ gg and gg ^ gg are excluded from the integral because gluons do not 
fragment into hadrons. The fragmentation functions used for this calculation are the 
original fragmentation functions of Feynman and Field [12] . For the pp — > tt X reactions, 
the fragmentation functions are 



Dl {z) = {z) = 
DT{z)=Dl\z) = 



l-z 



(^+l)(l_^)n^ 



(n + 1) (1 - ;.)", 



1 



D: {z)^Df{z)^D:{z)^(3' 



and for pp ^ K X reactions, the fragmentation functions are 



(n + l)(l -;.)", 
(n + l)(l-z)". 



D. 



K+ 



l-z 



(n + 1) {l-z)\ 



Dfiz) = Dfiz) = Dfiz) = Dfiz) = ^ 



Df {z)^Dr{z) 



D. 



K 



/3 + 



l-z 



, (z) = (z) = Df (z) = Df (z) = ^ 



^) {n + l){l-zr, 
n + l){l-zr, 

{n + l)(l-zr, 



where P — 0.4. The distributions of c, b and t quarks are sufficiently low that 

D'iiz) = D^{z) = D'^iz) = 0, 



(230) 

(231) 
(232) 
(233) 

(234) 

(235) 
(236) 

(237) 
(238) 



24 



for any hadron h. Feynman fixed n = 2 in his original paper. In this work, n is a 
parameter freely adjusted to fit data. There is a subtlety involved in summing all the 
parton contributions over a and b in Eq. (216) that is related to the relative probabilitistic 
nature of the parton distributions. The best way to explain it is by the way of an analogy. 
Suppose there is a room A occupied by 7 boys and 5 girls. The ratio of boys to girls in 
that room is 7:5, or : in normalized format. Suppose there is another room B of 14 
boys and 10 girls. The ratio is also 7:5. If a superintendent wants to know the number of 
children in room B, it is not enough to know the ratio of boys to girls. He must also be 
given a integral multiplicative constant. In this case, 2 times 7:5 or 24 times ^ : ^ gives 
the correct head counts which is 14:10. A similar situation arises in counting the number 
of partons in a hadron. The parton distributions are normalized to unity so that they give 
only the relative distributions of the partons. The parton distributions give only the ratios 
of the partons in a hadron but not their numbers. In order to sum over a and b partons 
in Eq. (216) correctly, one must be provided with an integral multiplicative constant for 
each of the hadrons A and B. These multiplicative integral constants cannot be known a 
priori but are determined a posteriori by fitting data. In other words, Eq. (216) can be 
modified as 

= ^^^^ y2 [ dxa I dxb 

^fA/a{Xa. Q') fB/biXb, Q') D^^{z,, Q') - ^, (239) 

ox 

where Na and Nb are the multiplicative constants corresponding to fA/a{xa,Q^) and 
fB/b{xb, Q^) respectively and Z^{a,6} is the sum over the parton types instead of a sum over 
the partons per se. li A = B, the overall multiplicative constant, N^Nb, is an integer 
square. If A ^ B, the overall constant, Na Nb, is still an integer. In the case of fitting 
the pQCD calculations to experimental data of the pp ^ n X reaction, a factor of 100 
is missing if one simply sums over the parton types. It implies that the multiplicative 
constant for the parton distributions of proton is Np — 10. The cross sections for K 
production is approximately half of that of tt production. It implies that the multiplicative 
constant for the pp ^ K X reaction is Np = 7. For the purpose parameterizing the shape 
of the kaon production cross section, the exact value of the scale is not required. Therefore 
Np in the pp — > KX reactions is arbitrarily set to be the same as that of pp — > nX reactions 
at Np — 10. This choice is adequate because fits to kaon experimental data are not being 
pursued due to a lack of experimental data. 

Figures 1 and 2 illustrate the features of E(P(j/dp^ plotted against px- Fig. 1 shows 
that the invariant cross section is visibly suppressed at high px- Fig. 2 shows that the 
suppression increases as 6cm moves farther away from 90° and that the suppression is 
symmetric around 9cm — 90°. Fig. 3 shows that the Monte Carlo code is plagued by nu- 
merical noise at excessively high ^/s. These illustrations summarize the global properties 
of the invariant production of cross sections of pions and kaons. 
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It is observed that the invariant production cross sections have the same basic shape 
regardless of the reactions, i.e. an exponentially decaying function of the form exp{—ax'^) 
at low pt and a suppression at high pt which drops off to zero before the edge of suppresion 
at Pt < The cross section is at its maximum at pt = and decreases monotonically 

in Pt- The Feynman-Field code used in this calculation assumes that = Ap^. In other 
words, by combining the previous two statements, the cross section is at its maximum 
at = and decreases monotonically in Q^. This observation indicates that hadron 
fragmentation is more favorable at low in that a parton preserves more kinetic energy 
to be made available for hadron fragmentation. It is also observed that the cross section 
is suppressed at high px and that the edge of suppression of the cross section is \/s/2 
at low Pt and gradually increases toward high pt- The reason for this phenomenon is 
mostly due to the choice of = ipT such that s > oi equivalently \/s/2 > pt- It 
turns out that the edge of suppression along pt is a function of y^. The comments made 
so far apply to any angle 9cm when ^/s is replaced by sin These observations will be 
incorporated into the parameterization described in the next section. 



3.5 Monte Carlo Integration 

A sample ANSI-C Monte Carlo program that calculates Ed^a/dp^ written in using the 
vegas routine is given below. The standard Monte Carlo integration packages are taken 
from Numerical Recipes [3]. They include vegas2.c, rebin.c and ran2.c. The Numeical 
Recipes packages need the header files nrutil.h, nr .h and nrutil. c to run. The CTEQ6 
package is downloadable from the web at http : / /www. cteq.org. The original package 
is written in f ortran but is translated to C using f2c. The translated CTEQ6 package is 
contained in the file Ctq6Pdf . c. The f 2c libraries must be installed in the system for the 
translator to work. When the code is compiled in unix/linux, the following line must be 
typed at the command prompt to link the f 2c libraries to the executable: 

$$ gcc -o pip pip.c -If 2c -Im 

where pip . c is the filename of the program. 

Below is a sample code that calculates the invariant production cross section of the 
pp TT^X reaction. 

#include <stdio.h> 
#include <math.h> 
#define NRANSI 

#include "/recipes/c_d/nrutil .h" 
#include "/recipes/c_d/nr .h" 
#include "/recipes/c_d/nrutil . c" 
#include "/recipes/c_d/vegas2 . c" 
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#include "/recipes/c_d/rebin. c" 
#include "/recipes/c_d/ran2.c" 
#include "Ctq6Pdf.c" 

#define pi 3.1415926535897932384626433 

#define lcLmbda2 0.16 /* part of alpha_s */ 

#define nf 5.0 /* number of flavors (topless) */ 

#define beta 0.4 /* part of the fragmentation function */ 

#define fac 3 . 893792923e-28 /* conversion from GeV-2 to cm''2 */ 

long idum; 

void vegasCdouble regn[] , int ndim, double (*fxn) (double [] , double), 
int init, unsigned long ncall, int itmx, int nprn, 
double *tgral, double *sd, double *chi2a) ; 

double fxn(double *x, double wt) ; 

doublereal ctq6pdf_ (integer *iparton, doublereal *x, doublereal *q) ; 
int setctq6_ (integer *iset) ; 
int readtbl_ (integer *nu) ; 
integer nextun_(void) ; 

int polint_ (doublereal *xa, doublereal *ya, integer *n, 

doublereal *x, doublereal *y, doublereal *dy) ; 
doublereal partonx6_ (integer *iprtn, doublereal *xx, doublereal *qq) ; 

double xl,x2,s,t,u,Q2,as2; 

float plabjd; 

integer *ipartona,*ipartonb; 
doublereal *q; 

MAIN__(void) 
{ 

int t 1 , t 2 , t3 , ndim=2 , i , imax , init=0 , ncall=5000 , itmx=5 ,nprn=0 ; 
int n=2*ndim; 
integer *iset; 

double pT,dpT,pTmin=l .0,xT,Th; 
float pTmax,tmp; 
double *regn; 
double alo,as,theta; 
double *tgral,*sd,*chi2a; 
FILE *fp; 

regn=dvector (1 ,n) ; 
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tgral=dvector (1 , 1) ; 
sd=dvector (1 , 1) ; 
chi2a=dvector(l, 1) ; 
ipartona=lvector (1 , 1) ; 
ipartonb=lvector (1 , 1) ; 
q=dvector (1 , 1) ; 

*iset=lL; 

f p=f open ( " out . dat " , " w" ) ; 

printfC'plab = "); 

scanf ("y„f ",&plab) ; 

printf ("theta = "); 

scanf ("y.f",&tmp); 

printf ("d = "); 

scanf ("y.f",&d) ; 

printf ("pTmax = "); 

scanf ("y.f",&pTmax); 

printf ("length (10 or 100): "); 

scanf ("y.d",&imax) ; 

tl=time(0) ; 



/* CTEQ6 for MSbar */ 



setctq6_(iset) ; /* initialization for CTEQ6 routines */ 

dpT= (pTmax-pTmin) / imax ; 
pT=pTmin ; 

theta=tmp*pi/ 180 . ; 
s=2 . 0*0 . 938272*plab ; 
Th=tan(theta/2.0) ; 

regn[3]=1.0; /* xa_ma:x */ 

regn[4]=1.0; /* xb_max */ 

t2=time(0) ; 

f or(i=0; i<=imax; i++) { 
Q2=4.0*pT*pT; 
q[l]=sqrt(Q2) ; 

alo=12 . 0*pi/ (33 . 0-2 . 0*nf ) /log (q2/lambda2) ; 

as=alo* (1 . 0- (306 . 0-38 . 0*nf ) /4 . 0/pi/ (33 . 0-2 . 0*nf ) *alo 

*log(log(Q2/lambda2))) ; 
as2=as*as ; 
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xT=2.0*pT/sqrt(s) ; 
xl=0.5*xT/Th; 
x2=0.5*xT*Th; 
t=-s*x2; 
u=-s*xl ; 

regn[l]=xl/(1.0-x2) ; /* xa_min */ 
regn[2]=x2/(1.0-xl) ; /* xb_min */ 
vegasCregn, ndim, fxn, init, ncall, itmx, nprn, 

&tgral[l], &sd[l], &chi2a[l]); 
f printf (f p , "%f \ty,e\n" ,pT , f ac*tgral [1] *100 . 0) ; 
pT+=cipT ; 

} 

t3=time(0) ; 

f close(fp) ; 

f ree_dvector (regn, 1 ,n) ; 
free_dvector(tgral,l,l) ; 
f ree_dvector (sd, 1,1); 
f ree_dvector (chi2a, 1,1); 
f ree_lvector (ipartona, 1,1); 
f ree_lvector (ipartonb, 1,1); 
f ree_dvector (q, 1,1); 
printf ("\ntime used for overheads 
printf ("\ntime used for integrals 

} 

double fxnCdouble *x, double wt) { 
double cs,z,Fbar,f ,tmp; 
double sh,th,uh; /* internal mandelstam vars */ 

double sh2,th2,uh2; 

double M1,M2,M3,M4,M5; /* inv amp */ 

double Du,Dd,Ds; /* fragmentation functions */ 

z=x2/x[2] +xl/x[l] ; 
sh=x[l] *x[2] *s; 
th=x[l]*t/z; 
uh=x [2] *u/z ; 
sh2=sh*sh; 
th2=th*th; 
uh2=uh*uh ; 



= '/od seconds\n",t2-tl) ; 
= y.d seconds\n\n\a",t3-t2) ; 
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Ml=4 . 0/9 . 0* (sh2+uh2) /tli2 ; 

M2=4 . 0/9 . 0* ( (sh2+uh2) /th2+ (sh2+th2) /uh2) -8 . 0/27 . 0*sh2/uh/th; 

M3=4 . 0/9 . 0* ( (sh2+uh2) /th2+ (th2+uh2) /sh2) -8 . 0/27 . 0*uh2/sh/th ; 

M4= (uh2+th2) /uh/th/6 . 0-3 . 0/8 . 0* (uh2+th2) /sh2 ; 

M5=-4 . 0/9 . 0* (uh2+sh2) /uh/sh+ (uh2+sh2) /th2 ; 

f = (d+1 . 0) *pow ( 1 . 0-z , d) ; 

Fbar=(1.0/z-1.0)*f ; 

Ds=Dd=beta*beta*Fbar ; 

Du=beta*f+Ds; 

if (x[2]<x[l]*x2/(x[l]-xl)) cs=0.0; 
else {. 
cs=0.0; 

/* su->s */ 
ipartona[l] =3L; 
ipartonb[l]=lL; 

cs+=ctq6pdf _(&ipartona[l] ,&x[l] ,&q[l]) 

*ctq6pdf_(&ipartonb[l] ,&x[2] ,&q[l] )*Ml*Ds; 

/* subar->s */ 
ipartona[l] =3L; 
ipartonb[l]=-lL; 

cs+=ctq6pdf_(&ipartona[l] ,&x[l] ,&q[l]) 

*ctq6pdf_(&ipartonb[l] ,&x[2] ,&q[l] )*Ml*Ds; 

/* us->u */ 
ipartona[l] =1L; 
ipartonb[l]=3L; 

cs+=ctq6pdf_(&ipartona[l] ,&x[l] ,&q[l]) 

*ctq6pdf_(&ipartonb[l] ,&;x[2] ,&q[l] ) *Ml*Du; 

/* usbar->u */ 
ipartona[l]=lL; 
ipartonb [1] =-3L ; 

cs+=ctq6pdf_(&ipartona[l] ,&x[l] ,&q[l]) 

*ctq6pdf_(&ipartonb[l] ,&x[2] ,&q[l] ) *Ml*Du; 

/ * ub->u */ 
ipartona[l] =1L; 
ipartonb [1]=5L; 
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cs+=ctq6pdf_(&ipartona[l] ,&x[l] ,&q[l]) 

*ctq6pdf_(&ipartonb[l] ,&x[2] ,&q[l] ) *Ml*Du; 

/* ubar->u */ 
ipartona[l] =1L; 
ipartonb [1] =-5L ; 

cs+=ctq6pdf _(&ipartona[l] ,&x[l] ,&q[l]) 

*ctq6pdf_(&ipartonb[l] ,&x[2] ,&q[l] )*Ml*Du; 

/* uc->u */ 
ipartona[l] =1L; 
ipartonb [1]=4L; 

cs+=ctq6pdf _(&ipartona[l] ,&x[l] ,&q[l]) 

*ctq6pdf_(&ipartonb[l] ,&x[2] ,&q[l] )*Ml*Du; 

/* ucbar->u */ 
ipartona [1] =1L ; 
ipartonb [1] =-4L ; 

cs+=ctq6pdf_(&ipartona[l] ,&x[l] ,&q[l]) 

*ctq6pdf_(&ipartonb[l] ,&x[2] ,&q[l] )*Ml*Du; 

/* sd->s */ 
ipartona [1]=3L; 
ipartonb [1]=2L; 

cs+=ctq6pdf_(&ipartona[l] ,&x[l] ,&q[l]) 

*ctq6pdf_(&ipartonb[l] ,&;x[2] ,&q[l] ) *Ml*Ds; 

/* sdbar->s */ 
ipartona [1] =3L ; 
ipartonb [1]=-2L; 

cs+=ctq6pdf_(&ipartona[l] ,&x[l] ,&q[l]) 

*ctq6pdf_(&ipartonb[l] ,&x[2] ,&q[l] )*Ml*Ds; 

/* sb->s */ 
ipartona [1] =3L; 
ipartonb [1] =5L ; 

cs+=ctq6pdf_(&ipartona[l] ,&x[l] ,&q[l]) 

*ctq6pdf_(&ipartonb[l] ,&x[2] ,&q[l] )*Ml*Ds; 

/* sbbar->s */ 
ipartona [1] =3L; 
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ipartonb [1] =-5L ; 

cs+=ctq6pdf_(&ipartona[l] ,&x[l] ,&q[l]) 

*ctq6pdf_(&ipartonb[l] ,&x[2] ,&q[l] )*Ml*Ds; 

/* sc->s */ 
ipartona [1] =3L ; 
ipartonb [1]=4L; 

cs+=ctq6pdf _(&ipartona[l] ,&;x[l] ,&q[l]) 

*ctq6pdf_(&ipartonb[l] ,&x[2] ,&q[l] )*Ml*Ds; 

/* scbar->s */ 
ipartona [1] =3L; 
ipartonb [1] =-4L ; 

cs+=ctq6pdf_(&ipartona[l] ,&x[l] ,&q[l]) 

*ctq6pdf_(&ipartonb[l] ,&x[2] ,&q[l] )*Ml*Ds; 

/* ds->d */ 
ipartona [1] =2L ; 
ipartonb [1]=3L; 

cs+=ctq6pdf_(&ipartona[l] ,&x[l] ,&q[l]) 

*ctq6pdf_(&ipartonb[l] ,&x[2] ,&q[l] ) *Ml*Dd; 

/* dsbar->d */ 
ipartona [1] =2L; 
ipartonb [1] =-3L ; 

cs+=ctq6pdf _(&ipartona[l] ,&x[l] ,&q[l]) 

*ctq6pdf_(&ipartonb[l] ,&x[2] ,&q[l] ) *Ml*Dd; 

/* db->d */ 
ipartona [1] =2L ; 
ipartonb [1] =5L ; 

cs+=ctq6pdf_(&ipartona[l] ,&x[l] ,&q[l]) 

*ctq6pdf_(&ipartonb[l] ,&x[2] ,&q[l] )*Ml*Dd; 

/* dbbar->d */ 
ipartona [1] =2L ; 
ipartonb [1] =-5L ; 

cs+=ctq6pdf _(&ipartona[l] ,&x[l] ,&q[l]) 

*ctq6pdf_(&ipartonb[l] ,&x[2] ,&q[l] )*Ml*Dd; 

/* dc->d */ 
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ipartona [1] =2L ; 
ipartonb[l]=4L; 

cs+=ctq6pdf _(&ipartona[l] ,&x[l] ,&q[l]) 

*ctq6pdf_(&ipartonb[l] ,&x[2] ,&q[l] )*Ml*Dd; 

/* dcbar->d */ 
ipartona [1] =2L; 
ipartonb [1] =-4L ; 

cs+=ctq6pdf_(&ipartona[l] ,&x[l] ,&q[l]) 

*ctq6pdf_(&ipartonb[l] ,&x[2] ,&q[l] ) *Ml*Dd; 

/* ud->u */ 
ipartona [1]=1L; 
ipartonb [1]=2L; 

cs+=ctq6pdf_(&ipartona[l] ,&x[l] ,&q[l]) 

*ctq6pdf_(&ipartonb[l] ,&;x[2] ,&q[l] ) *Ml*Du; 

/* udbar->u */ 
ipartona [1] =1L; 
ipartonb [1]=-2L; 

cs+=ctq6pdf _(&ipartona[l] ,&x[l] ,&q[l]) 

*ctq6pdf_(&ipartonb[l] ,&x[2] ,&q[l] ) *Ml*Du; 

/* du->d */ 
ipartona [1] =2L ; 
ipartonb [1]=1L; 

cs+=ctq6pdf_(&ipartona[l] ,&x[l] ,&q[l]) 

*ctq6pdf_(&ipartonb[l] ,&x[2] ,&q[l] )*Ml*Dd; 

/* dubar->d */ 
ipartona [1] =2L ; 
ipartonb [1]=-1L; 

cs+=ctq6pdf _(&ipartona[l] ,&x[l] ,&q[l]) 

*ctq6pdf_(&ipartonb[l] ,&x[2] ,&q[l] )*Ml*Dd; 

/* ■u.u->u */ 
ipartona [1] =1L; 
ipartonb [1]=1L; 

cs+=ctq6pdf _(&ipartona[l] ,&x[l] ,&q[l]) 

*ctq6pdf_(&ipartonb[l] ,&x[2] ,&q[l] )*M2*Du; 
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/* dd->d */ 
ipartona[l] =2L; 
ipartonb[l]=2L; 

cs+=ctq6pdf_(&ipartona[l] ,&x[l] ,&q[l]) 

*ctq6pdf_(&ipartonb[l] ,&x[2] ,&q[l] ) *M2*Dd; 

/* ss->s */ 
ipart ona [1] =3L ; 
ipartonb[l]=3L; 

cs+=ctq6pdf_(&ipartona[l] ,&x[l] ,&q[l]) 

*ctq6pdf_(&ipartonb[l] ,&;x[2] ,&q[l] )*M2*Ds; 

/* uubar->u */ 
ipartona[l] =1L; 
ipartonb[l]=-lL; 

cs+=ctq6pdf_(&ipartona[l] ,&x[l] ,&q[l]) 

*ctq6pdf_(&ipartonb[l] ,&x[2] ,&q[l] ) *M3*Du; 

/* ddbar->d */ 
ipartona [1] =2L ; 
ipartonb[l]=-2L; 

cs+=ctq6pdf_(&ipartona[l] ,&x[l] ,&q[l]) 

*ctq6pdf_(&ipartonb[l] ,&x[2] ,&q[l] )*M3*Dd; 

/* ssbar->s */ 
ipartona [1] =3L ; 
ipartonb [1] =-3L ; 

cs+=ctq6pdf _(&ipartona[l] ,&x[l] ,&;q[l]) 

*ctq6pdf_(&ipartonb[l] ,&x[2] ,&q[l] )*M3*Ds; 

/* gg->u,d,s */ 
ipartona [1] =0L; 
ipartonb [1]=0L; 

cs+=ctq6pdf _(&ipartona[l] ,&;x[l] ,&;q[l]) 

*ctq6pdf_(&ipartonb[l] ,&x[2] ,&q[l] )*M4* (Du+Dd+Ds) ; 

/* ug->u */ 
ipartona [1] =1L; 
ipartonb [1]=0L; 

cs+=ctq6pdf_(&ipartona[l] ,&x[l] ,&q[l]) 

*ctq6pdf_(&ipartonb[l] ,&x[2] ,&q[l] ) *M5*Du; 
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/* dg->d */ 

ipartona[l]=2L; 

ipartonb[l]=OL; 

cs+=ctq6pdf_(&ipartona[l] ,&x[l] ,&q[l]) 

*ctq6pdf_(&ipartonb[l] ,&;x[2] ,&q[l] ) *M5*Dd; 

/* sg->s */ 
ipartona[l] =3L; 
ipartonb[l]=OL; 

cs+=ctq6pdf _(&ipartona[l] ,&x[l] ,&q[l]) 

*ctq6pdf_(&ipartonb[l] ,&x[2] ,&q[l] )*M5*Ds; 

cs*=as2/sh2/z; 

} 

return cs; 

} 

#undef NRANSI 

4 Conclusion 

The Lund and Feynman-Field models are derived and implemented numerically in this 
work. New results based on these models will published as a separate paper. 

References 

[1] B. Andersson, The Lund Model, (Cambridge Univ. Press, 1998). 

[2] R. D. Field, Applications of Perturbative QCD, (Addison- Wesley, 1989). 

[3] W. H. Press, S. A. Tcukolsky, W. T. Vcttcrling and B. P. Flannery, Numerical 
Recipes in C, (Cambridge: Cambridge, 1997). 

[4] A. Tang, "Topics in Hadronic Physics," Ph. D. dissertation. Physics Department, 
University of Wisconsin-Milwaukee (December 2002). 

[5] J. Polchinski, Introduction to String Theory, (Cambridge Univ. Press, 1998). 

[6] B. Hatfield, Quantum Field Theory of Point Particles and Strings, (Addison- Wesley, 
1992). 



35 



[7] B. Andersson, S. Mohanty and F. Soderberg, hep-ph/0106185, (16 June 2001). 

[8] V. V. Abramov, et al, Nuclear Phys., B173, 348 (1980). 

[9] G. P. Lepage, Cornell University Publication, CLNS-80/447 (1980). 

[10] J. Pumplin, D. R. Stump, J. Huston, H. L. Lai, P. Nadolsky, and W. K. Tung, 

hep-ph/0201195 (22 January 2002). 

[11] S. R. Blattnig, S. R. Swaminathan, A. T. Kruger, M. Ngom and J. W. Norbury, 
NASA Technical Paper, NASA/TP-2000-210640, (December 2000). 

[12] R. D. Field and R. P. Feynman, Nuclear Phys., B136, 1-76 (1978). 

[13] B. L. Combridge, J. Kripfganz and J. Ranft, Phys. Lett, 70B, 2, 234-238 (26 
September 1977). 



36 



le-25 




4 5 6 
PT(GeV) 



10 



Figure 1 : A sample plot of the Feynman-Field model for pp — > t:^X at different energies 
and 0cm — 90°. The plot shows that the invariant cross section is suppressed at high pt- 
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Figure 2: A sample plot of the Feynman- Field model for pp — > tt'^X at different angles. 
The plot shows that the invariant cross section is symmetric around 9cm — 90° and is 
suppressed at 9cm 90°- 
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Figure 3: A sample plot of the Feynman-Field model for pp vr+X at excessively high 
energies and 9cm = 90°. The plot shows that the Monte Carlo code breaks down at 
ultra-high energies. 
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Figure 4: Comparisons of the Feynman-Field model fit for pp — ^ vr+X for various n at 
0cm = 90°. The references of the data sets Abromov 80, Abramov 81, Jaffe 89 and Busser 
76 are Nuclear Phys., B173 348 (1980), ZETFP 33 304 (1981), Phys. Rev., D40 2777 
(1989) and Nuclear Phys., B106 1 (1976) respectively. 
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Figure 5: The yoyo motion of a quark-antiquark pair confined by a linear potential. The 
two quarks are assumed to be massless. The mass square of the meson is proportional to 
the area of each square. The meson shown here is at rest. 
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Figure 6: The breakup of a quark- ant iquark pair along a surface of constant proper time 
r. The bold arrows represents the velocities of the produced mesons. The breakup points 
are labelled as vertices Vi to V^. 



42 




Figure 7: The geometry of the kinematics of two adjacent vertices Vi and V2. Vi and V2 
are vertices or spacetime positions which also represent the energy carried by the string 
field. The quark moves along the positive light-cone and the antiquark moves along the 
negative light-cone. The antiquark of Vi combines with the quark of V2 to produce a 
meson of mass m. W+i is the energy of Vi along x+. z+W+i is the fraction of energy 
used to create a quark from V2. W^2 is the energy of V2 along x_. z_W_2 is the fraction 
of energy used to create an antiquark from Vi. Ai, A2, and are the areas of the 
rectangles. The figure is taken from reference [1]. 
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Figure 8: The geometry of the kinematics of involving the total area of the spacetime 
diagram such that Afotai = T + Arest- The total energy square s — W+W^ is represented 
by the rectangle labelled s. p± is the energy of the parent quark (antiquark) along the 
± light-cone coordinates. ZQnP+ is the fraction of energy used to create a quark from Vn- 
w? / {zQnP+) is the energy used to create an antiquark from V^_i. The figure is taken from 
reference [1] with minor modifications. 
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Figure 9: Parton a of hadron A collides with parton b of hadron B producing partons c 
and d. The transverse momenta, pr, of c and d are equal and opposite. 
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